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We consider pure quantum states of N 3> 1 spins or qubits and study the average entanglement 
that can be localized between two separated spins by performing local measurements on the other 
individual spins. We show that all classical correlation functions provide lower bounds to this 
localizable entanglement, which follows from the observation that classical correlations can always be 
increased by doing appropriate local measurements on the other qubits. We analyze the localizable 
entanglement in familiar spin systems and illustrate the results on the hand of the Ising spin model, 
in which we observe characteristic features for a quantum phase transition such as a diverging 
entanglement length. 
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The mathematical description of multiparticle quan- 
tum systems plays an important role in several branches 
of physics. The main difficulty stems from the fact that 
the number of parameters needed to describe a quantum 
state grows exponentially with the number of particles. 
However, sometimes it is possible to capture the most 
relevant physical properties by describing these systems 
in terms of very few parameters. This is the case, for 
example, in quantum statistics, where two-particle cor- 
relations play a fundamental role. They allow us to un- 
derstand several complex physical phenomena, like phase 
transitions. Furthermore, they give rise to concepts like 
correlation length, which quantifies a very intuitive prop- 
erty of these systems. 

Multiparticle systems are also of central interest in 
the field of quantum information and, in particular, the 
quantification of the entanglement contained in quantum 
states. The reason is that entanglement is the physical 
resource to perform some of the most important quantum 
information tasks, like quantum information transfer (cfr. 
teleportation) or quantum computation. 

Given the common interest of quantum statistical me- 
chanics and quantum information in multiparticle sys- 
tems it is natural to try to describe the physical phenom- 
ena, like quantum phase transitions, appearing in (e.g.) 
spin systems from the point of view of entanglement. The 
main restriction one encounters is the fact that there ex- 
ist very few measures of multiparticle entanglement with 
a clear physical meaning. In any case, since entangle- 
ment measures correlations (for pure states), one would 
expect that reasonable entanglement measures be inti- 
mately connected to the correlation functions widely used 
in the context of quantum statistical mechanics. This is 
not the case, however, if one studies the behavior of the 
entanglement of formation between two separate spins 
after tracing out the rest Although this ap- 

proach exhibits a very pronounced (universal) behavior 
of this quantity at the transition point, it rapidly van- 
ishes as the distance between the spins goes beyond 3 or 
4, and thus it is not related to the correlations possessed 



by the state. Note also that the approach by Vidal et al. 
p, in which they study the scaling of entanglement of 
a block of spins with the other ones, does not quantify 
the entanglement as a resource (that could be used, for 
example, for teleportation). 

In this Letter we introduce a new concept which we call 
Localizable Entanglement (LE). On the one hand, this 
quantity has a very well defined physical meaning which 
treats entanglement as a truly physical resource. On the 
other, it establishes a very close connection between en- 
tanglement and correlation functions, as one would nat- 
urally expect. The LE of two particles is the maximal 
amount of entanglement that can be localized in these 
two particles, on average, by doing local measurements 
in the rest of the particles. The LE naturally leads to 
the definition of entanglement length, which measures the 
typical length scale at which the LE decays. The LE 
has an operational meaning which applies to situations 
in which out of some multiparticle entangled state one 
would like to concentrate as much entanglement as possi- 
ble in two particular particles. This occurs, for example, 
in the context of quantum repeaters |jj or in the con- 
text of quantum transport with spin systems (spintronics 
8J). The determination of the LE is a formidable task 
since it involves an optimization over all possible local 
measurement strategies, and thus cannot be determined 
in general. We have nevertheless managed to determine 
tight upper and lower bounds. The first ones stem from 
considering joint measurements on the rest of the spins, 
which is closely related to the concept of entanglement 
of assistance The second and more interesting ones 
can be derived by proving that there always exist local 
measurement which do not decrease the existing correla- 
tions between the two spins, something that despite its 
generality, up to our knowledge has not been considered 
in the context of quantum information. 

The usefulness of these findings will be illustrated on 
the hand of the entanglement present in the ground states 
of standard spin Hamiltonians. Let us however first con- 
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sider a simple example involving the A^-qubit GHZ-state: 
\GHZ) = ^=({00 ■■■0} + 



In this case, the LE is maximal as it is possible to create 
a Bell-state between two arbitrary qubits by measuring 
the other qubits in the |+), | — ) basis (here |±) = (|0) ± 
\l))/y/2). The existence of these quantum correlations 
could also have been revealed by studying the classical 
correlation functions 
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where i, j denote the positions of the spins under interest 
and a, (3 label the Pauli matrices. The correlation in the 
a = (3 = z-direction is the maximal possible one. It will 
indeed be shown that correlation functions yield lower 
bounds to the LE. However, in general the presence or ab- 
sence of classical correlations only gives a coarse-grained 
picture of the entanglement that ought to be created be- 
tween two distant spins. It is easy to find examples of 
highly entangled quantum states exhibiting no classical 
correlations whatsoever between any pair of spins. As 
an example, consider the so-called cluster states ^(J, ob- 
tained by the unitary evolution of an initially separable 
state under the action of the Ising Hamiltonian: 



^(^ 1 (|o)^ +1 ) + |i). 1 ))(|o) w 
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When N > 5, all reduced 2-qubit density operators are 
proportional to the identity and hence no correlations 
exist between any two spins. However, suitable local 
measurements on any N — 2 qubits can always create 
a Bell-state between the two remaining ones jlOj . hence 
indicating maximal LE. 

Let us next give a formal definition of the LE. Consider 
a pure state of N spins. Then the localizable entan- 
glement Eij(ip) is variationally defined as the maximal 
amount of entanglement that can be created (i.e. local- 
ized) , on average, between the spins i and j by performing 
local measurements on the other spins. More specifically, 
every measurement basis specifies a pure state ensem- 
ble £ — {p s , \(/> s )} consisting of at least 2^ N ~ 2 ^ elements 
counted by the index s. In this notation p s denotes the 
probability to obtain the two-spin state \<p s ) after per- 
forming the measurement |s) on the assisting spins of 
our N partite spin system. The LE is then defined as 



Eij = max) p s E{\<f> a )), 



the maximal correlation function 0] for a pure state 
of two qubits coincides with the concurrence, given by 
C(\<ip) = a|00) + 6|01) + c|10)+d|ll)) = 2\ad-bc\. Note 
that the behavior of the LE as defined in terms of the 
entropy of entanglement would be very similar [l9| . 

Due to the variational definition, the LE is very diffi- 
cult to calculate in general. Moreover, in typically large 
spin systems one does not have an explicit parametriza- 
tion of the state under interest, but just information 
about the classical 1- and 2-particle correlation functions 
(which parameterize completely the 2-qubit reduced den- 
sity operator). It would therefore be interesting to derive 
tight upper and lower bounds to the LE solely based on 
this information. The upper bound can readily be ob- 
tained using the concept of entanglement of assistance 
9] , which would correspond to the LE if global or joint 
measurements were allowed on the other spins. It was 
shown in 0] that the entanglement of assistance can be 
calculated as follows: given the reduced density operator 
Pij and a square root X, pij = XX', then the entan- 
glement of assistance Eij(tp) as measured by the concur- 
rence is equal to the trace norm Tr\X T (a y ® <r y )X\. 

The lower bound is more subtle, and will be shown 
to follow from the following interesting theorem: given a 
(pure or mixed) state of N qubits with classical correla- 
tion Q l ^a between the spins i and j and directions a, f3, 
then there always exist directions in which one can mea- 
sure the other spins such that this correlation do not de- 
crease, on average. This automatically implies that there 
always exist local measurements that increase (or keep) 
the classical correlations. Surprisingly, this very general 
theorem seems not to have been noticed before, and is 
interesting on its own. It could be very useful in the 
context of cryptography (e.g. multipartite distribution 
of common randomness |lfil|L 

Let us next proof this theorem. Note that it is suf- 
ficient to consider mixed states of three qubits. Let us 
parameterize a mixed 3-qubit density operator by four 
4x4 blocks 

r Pi o- 

ct^ p 2 

Without loss of generality, let us consider the Q^j, cor- 
relations. The original correlations are completely de- 
termined by the diagonal elements of the density op- 
erator pi + p%. A von Neumann measurement in the 
|±) := cos(0/2)|O) ± sin(0/2)exp(±i0)|l) basis on the 
third qubit results in the hermitian unnormalized 2-qubit 
operators 



where E(\<f> s )) denotes the entanglement of \<f> s ). As 
we deal with pure states of two qubits, all entangle- 
ment measures are essentially equivalent, and in order 
to make the connection with correlation functions, we 
will measure the entanglement on the hand of the con- 
currence 01 : indeed, it can readily be checked that 



„ Pi + P2 , fa s Pl - P2 
X± = — - — ±COS(0) 



i(a — a') 



±sin(0) ( cos(0) — h sm(<^>) ) 



Defining p± = Tr(X±), we have to prove that there al- 
ways exist parameters 0, cf> such that the correlation, on 
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average, does not decrease: 

p+ \q zz (x + Ip+)\+p-\QUX-/p-)\ > \Q„(X+ + X-)\. 

Without loss of generality we can assume that a — 
Q ZZ {X + + X_) is positive. Removing the absolute value 
signs and parameterizing 

x := [cos(0); sin(0) cos(</i); sin(#) sin(</>)], 

some straightforward algebra yields the sufficient inequal- 
ity 
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where c is such that p± = (1 ± c T x)/2, and f3, Q are 
defined as 3 x f and 3x3 blocks of the matrix 



S = R T (a y a v )R = 



a p T 

P Q 



Here R is the real 4x4 matrix whose columns consist of 
the diagonal elements of the matrices (pi + P2), (pi — P2), 
(o- + o-t), i(a-at). 

Due to Sylvester's law of inertia [12j, we know that 
S has two positive and two negative eigenvalues. Now 
Q — P(3 T /a is the Schur complement of a, and hence cor- 
responds to a principal 3x3 block of the matrix S . 
Due to the interlacing properties of eigenvalues of prin- 
cipal blocks 0] > it follows that it has either two positive 
and one negative eigenvalue or two negative and one pos- 
itive one. And this of course ensures that there always 
exists an x such that the inequality Q is fulfilled, com- 
pleting the proof. □ 

Note that the proof is constructive and allows to deter- 
mine a measurement strategy that would at least achieve 
the bound reported. Note also that exactly the same 
proof applies when correlations between larger blocks of 
spins would be considered [l3| . 

Let us now show how this theorem yields a lower bound 
to the LE. Given an initial pure state of N qubits, we 
know that the measurement of the first, second, ... N — 
3'th qubit can be chosen such that on average the final 
correlations do not decrease. But we end up with a pure 
state of two qubits, for which the concurrence is equal 
to the maximal correlation. A lower bound to the LE 
is therefore given by the maximal correlation function, 
and following the previous theorem, there is always a 
constructive way of determining a measurement strategy 
that achieves this lower bound. Surprisingly, we will see 
that this lower bound seems to be the exact value for the 
LE in the case of many systems of interest. 

The previous findings can readily be applied to the 
study of entanglement in translational invariant ground 
and excited states of spins arranged in a regular lattice. 
In quantum statistics and more specifically in the study 



of quantum phase transitions, the correlation length is 
the canonical parameter of interest. The concept of LE 
readily lends itself to define the related entanglement 
length £e as the typical length scale at which it is possible 
to create Bell states by doing local measurements on the 
other spins. More specifically, the entanglement length 
is finite if and only if the the LE Ei y i +n ~ exp(— u/^e) 
decays exponentially in n, and the entanglement length 
£s is defined as the constant in the exponent in the limit 
of an infinite system (see also Aharonov 
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Note that a diverging correlation length automatically 
implies a diverging entanglement length and hence long- 
range quantum correlations, although the converse is not 
necessarily true|20j. 

More specifically, we have studied ground states of 
two spin interaction Hamiltonians with parity symmetry 
®iLi& l z of the form: 



T J (7 G 
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with coupling coefficients 7^ > 7*-' > (note that this 
model also includes spin systems in higher dimensional 
lattices) . Extensive numerical calculations on systems of 
up to 20 qubits showed that our lower bound is always 
very close to the LE, and typically is even equal to the 
exact value of the LE 0|: this is surprising and high- 
lights the power of the given lower bound. Note also 
that whenever the parity symmetry is present, the upper 
and lower bounds can easily be calculated in terms of the 
correlation functions: 



maX (Q % x X , Qyy , 
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As an illustration, we consider the ground state of the 
Ising Hamiltonian (7^ = XS atX Sj t i + i; 7* = 1), which has 
been solved exactly |l7| and exhibits a quantum phase 
transition at A = 1. In this case, the maximal classical 
correlation is always given by Q xx , and we conjecture 
that it is equal to the LE (for more details, we refer to 
(l^L For A < 1, the LE is small as the ground state is 
almost separable, and the entanglement length is finite 
(see Fig. 1) as the LE decreases exponentially with the 
spin distance. At the quantum critical point A = 1, the 
behavior of the LE changes drastically, as it decreases 
polynomially Ei^ +n ~ rt -1 / 4 , thus leading to a diverging 
entanglement length In the case A > 1 we also get 
£e = 00, as the LE saturates to a finite value given by 
M x = 1/4(1 — A -2 ) 1 / 4 . Indeed, the ground state is then 
close to the GHZ-state. This behavior is illustrated in 
Fig. 1. 
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FIG. 1: Localizable Entanglement Eij and Correlation func- 
tion Q^ x (solid) and upper bound (dashed) as a function of 
the coupling parameter A for the infinite Ising chain with spin 
distance n = 1 (left) and n — 10 (right). 




FIG. 2: Correlation function Q^ +4 (solid), Localizable En- 
tanglement Ei.i+4 (diamonds), and its upper bound (dashed) 
as a function of the coupling parameter A for the ground state 
of the finite Ising chain (N = 14) in the case of broken par- 
ity symmetry by a small perturbating magnetic field in the 
x-direction. 



In a more realistic setup however, the parity symmetry 
of the Ising Hamiltonian will be broken by a perturba- 
tion and the ground state for large coupling will also be 
separable. Indeed, the energy gap between the lowest en- 
ergy states with different parity decays exponentially in 
the number of qubits in the region A > 1 , and henceforth 
the ground state becomes a superposition of these two 
states. The calculation of its LE and the corresponding 
lower and upper bounds are depicted in Fig. 2. 

These results give a clear illustration of the intimate 
connection between classical correlations and entangle- 
ment in the case of ground states of translational invari- 
ant Hamiltonians. The plethora of results concerning 
classical correlation functions, such as diverging correla- 
tion length at quantum phase transitions, can now be 
interpreted from the perspective of quantum information 
theory; one could argue that the status of classical cor- 
relations has been lifted to the one of (useful) quantum 
correlations. 

In conclusion, we have introduced the notion of local- 
izable entanglement. It has a nice operational meaning 
as it quantifies the amount of useful entanglement that 
can be created between two spins by doing local mea- 
surements on all other spins. We proved that classical 



correlation functions always provide lower bounds to the 
LE, showing that the presence of classical correlations 
is sufficient to be able to create Bell-like quantum cor- 
relations. As a side-product, we proved that classical 
correlations in multipartite mixed quantum states can 
always increase by doing appropriate measurements on 
the other qubits. Finally, we demonstrated the useful- 
ness of these concepts in the context of spin systems on 
a lattice, provided a natural definition of entanglement 
length, and showed that it diverges at a quantum phase 
transition. Further generalizations and applications will 
be presented elsewhere [l3|. 
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The basis with maximal possible correlation can in gen- 
eral be found by calculating the singular value decompo- 
sition of the 3x3 matrix Q a p, and the largest singular 
value corresponds to the maximal possible correlation. 
For a given concurrence C, the entropy of entanglement 
is given by /(C) = H ((1 + vT - C 2 )/2) with H(x) the 
Shannon entropy. Due to the convexity of /(C), the 
following bounds hold: /(I] s PsC s ) < ]C s p s /(C s ) < 

This leaves open the possibility of Hamiltonians exhibit- 
ing a phase transition as indicated by a diverging entan- 
glement length at the critical point but with a correlation 
length remaining finite. 



